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Abstract. The — 15.1 years period found in the long-term UBV photoelectric and photographic photometry of 
the symbiotic nova V1016 Cyg is detected also in the (J — K) colour index and in the UV continuum and emission 
line fluxes from IUE and HUT spectra. It could be interpreted either as the effect of recurrent enhanced mass loss 
episodes from the Mira type variable companion to a hot component along its ultra-wide orbit (proposed from 
recent HST observations) or the true orbital period of the inner, unresolved binary of a triple system. A 410-day 
delay of the maximum of UV emission lines fluxes with respect to the maximum of continuum was found. The 
pulsation period of the Mira type variable was improved to 474±6 days. 
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1. Introduction 

V1016 Cyg (MHa 328-116) is a member of a small sub- 
group of symbiotic stars, called symbiotic novae, also in- 
cluding V1329 Cyg and HM Sge, whose outburs ts lea d 
to a nebular spectrum (Miirset & Nussbaumer 1994 ). 
Symbiotic novae are wide interacting binaries, where mat- 
ter from a late-type giant is transferred onto the surface 
of the more compact companion. The nova-like optical 
outburst (Am — 5 — 7 mag), lasting decades, is caused 
by a thermonuclear runaway on the surface of a wind- 
accreting white dwarf after the critical amount of material 
has been accumulated (cf. Mikolajewska & Kenyon, 1992Q . 
V1016 Cyg underwent such nova-like outburst in 1964 
(McCuskey |1965 ) . The object is classified as a D-type sym- 
biotic, the cool component being a Mira type variable em- 
bedded in a dust envelope whose pulsation period turned 
out to be -478 (Munari |1988|). The onset of a dust for- 



mation episode in 1983 is reported by Taranova & Yudin 

( Eggp - 

The orbital perio d of V 1016 Cyg is not yet established. 
Taranova & Yudin (1983) made use of the increase of 
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* Partly based on observations obtained with the 
International Ultraviolet Explorer (IUE) satellite retrieved 
from the IUE Newly Extracted Spectra (INES) Archive and 
the Hopkins Ultraviolet Telescope (HUT) retrieved from the 
Multimission Archive at STScI (MAST). 



Balmer emission lines in combination with the appearance 
and disappearance of Fell lines to derive an orbital period 
of —20 years. Afterwards, Nussbaumer & Schmid ( 1988 ), 
though unable of recording two consecutive maxima, pro- 
posed an orbital period of 9.5 years on the basis of the 
apparent periodicity seen in the flux of OI and Mgll UV 
emission lines by means of the IUE satellite. At the same 
time Munari (1988), by resorting to IR observations taken 
over two decades, proposed instead a 6-year orbital period 
by modeling the sequence of dust obscuration episodes, 
likely related to the passage of the Mira type variable at 
the inferior conjunction in the system. 

Much longer periods have been proposed by 
Wallerstein fll98$ and Schild & Schmid fll996j ). In the 
former paper the author, assuming that the sharp Fell 
emission lines are formed in the chromosphere of the cool 
star so as to reflect its orbital motion, concludes that their 
observed radial velocities between 1978 and 1985 limit any 
high inclination orbit to a period greater than 25 years or 
to a large eccentricity. The latter analysis, based on spec- 
tropolarimetric data taken from 1991 to 1994, indicates 
that the orbital period is about 80 ± 25 years, though 
later observations obtained in 1997 put this result into 
question (Schmid, |1998[ ). Finally, Brocksopp et al. fl2002| ), 
adopting a projected separation of the two stellar compo- 
nents as large as 84 AU on the basis of their HST/WFPC2 
images, are forced to propose an astonishingly long orbital 
period of —544 years. 
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Fig. 1. Photographic (left) and UBV (right) light curves of V1016 Cyg. Arrows a,b,c and d mark the epochs of 
subsequent activity episodes of the system (see text). 



In the context of this still-to-be-settled issue, one has 
also to interpret the evidence for the 15-year periodic ac- 
tivity (not necessarily representing the soug ht orb ital pe- 
riod) recently presented by Parimucha et al. ( 200C| ) (here- 
after PAC) who gathered and analysed long-term pho- 
tographic, photoelectric and visual photometry of the ob- 
ject. The aim of the present paper is indeed to give further 
(i.e., multiwavelength) support to the existence of such a 
periodicity making use of both IR photometry and UV 
spectroscopy, as well as to investigate its own origin. 



2. Long-term photometry 

2.1. UBV photoelectric and photographic data 

The historical light curve of V1016 Cyg based on the pho- 
tographic and UBV data given in PAC is presented in 
Fig. 1. The light curve suggests four stages of activity 
marked by arrows in the figure: the pre-outburst flare a 
in 1949, the main nova-like outburst b in 1964 and two 
post-outburst, decreasing- amplitude flares c and d in 1980 
and 1994, respectively. Evidently the activity episodes af- 
fecting the system repeat themselves (though at quite a 
different intensity level) with an interval of «15 years. The 
ephemeris for the activity maxima calculated in PAC is as 
follows: 



JD-p^ = 2427590 
±250 



5510 xE. 
±90 



(1) 



It is worth noticing that both maxima recorded in 1980 
and 1994 followed appreciable brightness decreases which, 
in turn, could be interpreted as signatures of an enhanced 
mass transfer from the cool to the hot component. A sim- 
ilar effect has been detected in the light curve of t he ver y 
slow classical nova V723 Cas (Chochol & Pribulla, 1998[) . 



2.2. JHK data 



Infrared photometry of V1016 Cyg was published by 
Harve y (|1974| ), Kcnyon & Gallag her (|1983| ), Lor cnzett i 
et al. ( |1985D , Tara nova fc Yudin ( |1986| ), Mim ari (|1988D , 
Ananth & Leahy ( [1993D Kam ath & Ashok ( [l999h and 
Taranova & Schcnavrin (200C). We make use here of this 
whole dataset plus a few unpublished JHK observations 
obtained in 1993 (kindly provided by F. Strafella) to newly 
estimate the period of pulsation of the Mira type variable. 
The Fourier period analysis (Fig. 2) applied to the JHK 
data leads to the values 476±10, 472±11 and 474±10 days 
for J, H and K photometry, respectively. The adopted 
period of the Mira type variable pulsations (whose corre- 
sponding phase diagrams are also presented in Fig. 2) is 
thus P = 474±6 days. 

The long-term behaviour of the (J — K) color index of 
V1016 Cyg using all available infrared data is presented in 
Fig. 3. As pointed out by Whitelock fll987|) , the (J - K) 
color index in symbiotic Miras is little affected by the Mira 
type variable pulsation, being conversely very sensitive to 
the amount of circumstellar dust around the cool com- 
ponent. As a consequence, the abrupt change (~1 mag) 
of the IR color recorded in 1983 by Taranova & Yudin 
( 1986| ) comes likely from a short (though strong) dust for- 
mation episode. Although orbitally-related dust obscura- 
tion episodes are indeed expected in these systems, we 
interpret this unique, major episode which occurred three 
years after the observed maxima both in the U passband 
and space-borne UV (see below) as dust formation in the 
ejecta of the symbiotic nova (cf. Bode 1995). 



According to this view, one has to be cautious when 
interpreting possible periodicities emerging from the anal- 
ysis of the long-term behaviour of the (J — K) color index 
(Fig. 3) which shows an additional maximum in 1988 and 
a clear minimum in 1992, together with a more uncer- 
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Fig. 2. Upper panel: the Fourier power spectra of the pul- 
sation period of the Mira type variable in the JHK data. 
Lower panel: the resulting light curves phased with the 
ephemeris iT) max = 2447442 + 474 x£. The AH res and 
AK res are shifted by 1.5 and 3.0 mag, respectively. 



tain maximum (in 1973) and a further possible minimum 
(in 1977). Indeed, the phenomena lying beneath the two 
possible periods (namely, 2070 and 5630 days; see Fig. 4) 
one can identify by means of the Fourier analysis of its 
residuals (after removing a parabolic trend) is likely quite 
different. While the longer period, besides being close to 
the UBV photometric estimate, could well reflect the real 
orbital motion of the system, as suggested by the wave- 
like variation exhibited by the (J — K) residuals phased 
with the ephemeris (1) (see Fig. 4), the shorter period of 
5.6 years, close to the 6-year period interpreted by Munari 
( |l988|) as the orbital period of the binary system, is likely 
due to the superposition of several effects such as dust 
formation induced by the activity of the hot component 
and possible mass tranfer events from the cool to the hot 
component. 

3. UV spectroscopy 

The UV dataset discussed here includes both low and high 
dispersion IUE (International Ultraviolet Explorer) spec- 
tra spanning the interval from Aug. 1978 to Aug. 1995 
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Fig. 3. (J - K) color index of V1016 Cyg 



(for a complete list see Parimucha ( 2002| ), and a sin- 
gle HUT (Hopkins Ultraviolet Telescope) spectrum taken 
on March 6, 1995. Properly calibrated IUE NEWSIPS 
data have been retrieved from the IUE Newly Extracted 
Spectra (INES) Archive, while the calibrated HUT spec- 
trum was extracted from the Multimission Archive at 
STScI (MAST). 

Average continuum fluxes in emission-free, 20 A bins 
centered at A 1840 A and A 2530 A have been derived 
from low dispersion, large aperture IUE spectra (SWP and 
LWR/LWP cameras, respectively), plus the HUT spec- 
trum (for the A 1840 A region alone). Both measured and 
dereddened fluxes are given in Table 1 (having adopted 
E(B — V) — 0.28, according to Nussbaumer & Schild 
1981). Dereddened continuum fluxes are also shown in 
Fig. 5. 

High dispersion, large aperture IUE/SWP spectra 
were used to derive fluxes of OIII] (AA 1661,1666 A), NIII] 
(A 1750 A), Silll] (A 1892 A) and CIII] (A 1909 A) emis- 
sion lines. A single measurement of the OIII] lines flux 
was obtained also from the HUT spectrum of 1995. The 
emission line profiles show a broad component superposed 
to the dominant, narrow component (cf. Kindl, Marxer 
& Nussbaumer 1982). Fluxes given here are the result 
of fitting Gaussian profiles to the latter narrow compo- 
nent. Dereddened fluxes for OIII], NIII] and Silll] lines 
are shown in Fig. 6 and listed (together with the original 
fluxes) in Table 2. 

For the strong CIII] A 1909 A which falls by chance in 
the overlapping region between the IUE short-wavelength 
(SWP) and long-wavelength (LWR/LWP) camera spec- 
tral ranges, simultaneous flux estimates obtained from 
both cameras are shown in Fig. 7 and listed in Table 3. In 
order to make the best use of the whole dataset, Fig. 7 in- 
cludes also — as a lower limit — the estimated flux of CIII] 
in case the line saturation was marginal (i.e. affecting no 
more than 2 high-resolution pixels). 

It should be stressed that the availability of a large 
aperture spectrum taken on August 31, 1978 (namely, 
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Table 2. Measured (I) and corrected (II) (with E(B - V)=0.28) UV emission line fluxes in 1CT 12 erg cm" 2 s" 1 
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Note: a - HUT spectrum 



LWR 2228), allowing the measurement of the CIII] A 
1909 A emission line flux, does assure that the flux rise 
exhibited by the UV continua from mid-1978 to mid-1979 
(see Fig. 5) also affects the UV emission lines so as both 
UV continua and emission lines reflect the 1980 peak of 
activity of the system recorded in other wavebands (see 
above) . Strictly speaking, this verification turns out to be 
not possible for the remaining set of UV emission lines, 
as high resolution SWP spectra assured earlier than June 
24, 1979 were obtained through the non-photometric IUE 
small aperture. 

One can easily recognize that the evolution of both UV 
continua and emission lines from 1978 to 1996 is charac- 
terized by two maxima matching the activity epochs c,d 
shown by optical photometry. Moreover, interesting clues 
to the physical properties of the system could come from 
the observed delay of the onset of emission line activity, 
whose first maximum (at epoch JD 2444635 ±15) appears 
shifted by 410 days from the corresponding epoch for both 
UV continua (JD 2444225 ± 12). 

Finally, it should hold reader's interest the fact that 
the major dust formation episode discussed in the Sect. 2.2 
reflects in a recognizable flux drop of the majority of UV 
emission lines in 1983-84 (see Fig. 6). 



4. Discussion 

No doubt the symbiotic nova V1016 Cyg includes a pul- 
sating Mira type variable and an accreting white dwarf 
that underwent a thermonuclear outburst in 1964, leading 
straight to a nebular spectrum (FitzGerald et al. 1966 ). 
According to the ionization model of symbiotic binaries, 
the hot luminous component ionizes the neutral wind of 
the giant giving rise to the nebula in the system. Such a 
mechanism is confirmed by the strict coincidence of the 
epochs of maxima in the U passband (observed in 1980 
and 1994) and those recorded in the space-UV continuum 
combined with the 410-day delay of the OIII], CIII], NIII] 
and Silll] emission line maxima which, in turn, assures 
that the lines arc collisionally excited in the surrounding 
nebula when the fast wind from the hot object interacts 
with the slow wind from the Mira type variable. 

At present, the above well- understood phenomenology 
leaves still open two alternative scenarios for V1016 Cygni. 

If this symbiotic nova is a wide-orbit binary with a 



544-year period as proposed by Brocksopp et al. (2002), 
one is forced to interpret the periodic (15-yr) variations 
of the optical and UV continuum as induced by flares of 
the hot component triggered by the recurrent enhanced 
mass loss episodes from the M ira ty pe variable companion. 
According to Fleischer et al. ( |l995| ), this mass loss is most 
pronounced every few periods of the Mira type variable 
pulsations. It can trigger individual flares in the accretion 
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Table 1. Measured (I) and corrected (II) (with E(B 
V)—0.28) UV continuum flux in 10~ 14 erg cm s _1 A -1 
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Fig. 4. Upper panel: the Fourier power spectrum for the 
(J—K) residuals in the interval 1000 to 10000 days. Lower 
panel: the phase diagram of the above data corresponding 
to the ephemeris (1). A trigonometric polynomial has been 
adopted for the fit. 



disk of the hot white dwarf. The existence of the disk is 
supported by the 3D simulation of the wind accretion by 



the compact star (Theuns & Jorissen, 1993). 

Alternatively, V1016 Cygni could be interpreted as a 
triple system, whose inner, unresolved binary has an or- 
bital period of 15 years. If this is the case, flares would be 
triggered by an enhanced mass transfer from a cool giant 
during its periastron passage on the 15-year orbit. In this 
respect, a hint of enhanced mass transfer, coming just be- 
fore the weak flare in 1994, can be found in the long-term 
infrared photometry performed by Taranova & Schenavrin 
(2000) showing a brightness increase in the J and H pass- 
bands in 1992. Moreover, the coincidence of the observed 
maximum of the J — K index in 1988 (Fig. 3) and the 
wide minima of the UV continua (Fig. 5) would constrain 
to that epoch the inferior conjuction of a cool giant along 
its 15-year orbit. 
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Note: a - HUT spectrum 
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Fig. 6. Dereddened UV emission line fluxes 
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Table 3. Measured (I) and corrected (II) (with 
E(B - V)=0.28) CIII]1909A emission line fluxes in 
10 -11 erg cm -2 s _1 from SWP and LWP/R spectra 
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Note: 1 - lower limit 



